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Abstract 

The two spin-channel model is generalized to the case of transport of ferromagnetic excitations 
in electric conductors and insulators. The two channels are defined by reducing the ferromagnetic 
degrees of freedom to a bivaluated variable, i.e. to an effective spin one-half. The reduction is 
performed after defining the local magnetic configuration space by a sphere T, x , and integrating 
the relevant physical quantities over the two hemispheres £i and ££. The configuration space 
is then extended to the x direction for non-uniform magnetization excitations. The transport 
equations for both magnetic moments and magnetic energy are deduced, including the relaxation 
from one channel to the other. The heat transport equations for ferromagnets is deduced. 
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FIG. 1: Illustration of the Ti x configuration space of ferromagnetic degrees of freedom at point 
x. The separation in two hemispheres T, x = x£ U T, x allows the two effective spin-channels to be 
defined. The arrow in the center represents a mean value of the magnetization in the configuration 
space. 

In the context of spintronics, the effect of spin-injection and spin-accumulation is easy to 
describe on the basis of the two spin-channel model for electric carriers, if the magnetization 
is locally defined by the microscopic spin one-half s = ±h/2 [TH5]. In that case, the electric 
carriers define without ambiguity the two channels at any points x of the material: one 
channel for the up spin h/2 and the other channel for the down spin —h/2, with a fixed 
quantification axis. However, in the case of ferromagnetic interactions between electric car- 
riers (e.g. for spin- waves excitations), this definition is a-priori not valid since the magnetic 
carriers are not longer defined locally by a spin one -half. Beyond, in the case of electric 
insulators, the definition of two spin-channel seems to be problematic because it cannot be 
based on an assembly of delocalized quasi-particles. 

We show in this paper that the two spin-channel model can nevertheless be generalized to 
any kind of macroscopic ferromagnetic excitations (quantal or classical). This generalization 
is based on a reduction method that allows the continuous magnetic degrees of freedom to 
be reduced to a bivaluated variable at each point x of the usual space, i.e. to a local effective 
"spin" one-half. The reduction is performed after defining the local magnetic configuration 
space over the sphere S x (section I below), and integrating the relevant physical quantities 
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over each hemisphere (see Fig. 1) for each point x of the usual space (section II). The config- 
uration space is then extended to the x direction for non-uniform magnetization excitations, 
and the transport equations are deduced (section III). 

This generalization of the two spin-channel model shows that the effects generated by 
spin-polarized heat currents are similar to that generated by electric spin-polarized currents. 
Typically, the spin-accumulation effect occurring in spintronics devices can be generalized to 
thermal spin- accumulation [7]. The similarity of the transport equations of heat - together 
with the usual Seebeck and Nernst effects - could explain the so-called spin-Seebeck and 
spin-Nernst effects observed recently on various materials [HHTi]. 

A. Definition of the localized ferromagnetic configuration space T, x 

In this section, we focus on the localized uniform ferromagnetic moment M(x) = M s e r (x) 
defined with radial unit vector e r (x) and magnetization at saturation M s at point x. The 
approach used is the mesoscopic non-equilibrium thermodynamic theory (MNET) [T5l [16] 
applied to rotational brownian motion [T7J [IB] . 

In order to treat statistically the ferromagnetic degrees of freedom, a statistical ensemble 
of a large number of ferromagnetic moments is defined on the configuration space H x (Fig. 
1). Each magnetic moment is described by its position {9,(p} on the sphere H x of radius 
M s . The angle 9 is associated to the radial unit vector eg and the angle (p is associated to 
the azimuth unit vector e v . 

The statistical distribution of the magnetic moments on the sphere S x is then defined by 
the density p F (9,<p;x) per units of solid angle dfl = sin9d9d<p [19] and per unit length 5x. 
The function p F (9,ip) is a solution of the rotational Fokker-Planck equation derived at the 
end of this subsection. At equilibrium, the density is given by the Boltzmann distribution 
defined by the ferromagnetic potential V F (9,(p;x), and such that the effective magnetic 
field vanishes H e ff = — W F = (H e ff contains all deterministic contributions: external 
magnetic field, dipolar field, anisotropy field, exchange field, etc). For out-of-equilibrum 
states, fluctuations and diffusion plays a fundamental role. The distribution is given by the 
ferromagnetic chemical potential that takes the form fi F = kTln{p F ) + V F [TBI [T9l 120] . 
The generalized force H e ff + h = -Vs/i F produces a current of ferromagnetic moments 
J F (9,(p;x) = p F du r /dt, that is flowing on the surface of the sphere. This generalized force 
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contains a deterministic part H e ff = — W F and a diffusive part h = {kT / p F )V p F \TT\ 120]. 
If the ferromagnetic system is closed, the magnetic moments are conserved on E x so that: 



dp 



F 



div^J F . (1) 



dt 

On the other hand, the ferromagnetic energy u F is also described by a the continuity 
equation. However, the system under interest is not adiabatic and the conservation equation 
for the ferromagnetic energy takes the form: 

du F 

= -diwJ u + r (2) 

where r describes the dissipation in the environment and J u is the flux of ferromagnetic 
energy in the configuration space H x . If the temperature T(x) is uniform over the config- 
uration space S x (i.e. at each point x), the flux of heat is zero, and the power dissipated 
by the ferromagnetic system is reduced to the effect of the flux of magnetic moments: the 
ferromagnetic power dissipated or stored at each point x of the ferromagnet reads 

v F = - j P.v^dn (3) 

The application of the second law of thermodynamics allows the transport equation to be 
deduced by writing the relation that links the generalized flux to the generalized force. 
Both quantities, flux and forces, are then related by the Onsager matrix of the transport 
coefficients I 

J F = -CV^p F (4) 

where the flow J F is a two component vector defined with the unit vectors {e v ,eg} of 
H x . Accordingly, the Onsager matrix is a 2x2 matrix defined by four transport coefficients 
{Lgg, L0 V , L^o, L w }. The Onsager reciprocity relations impose that Lg v = —L v g. Further- 
more, assuming that the dissipation is isotropic, we have L ee = L w . Introducing a as the 
ratio of the off-diagonal to the diagonal coefficients; a = Lg^/Lgg, the ferromagnetic kinetic 
equation is defined by two ferromagnetic transport coefficients L F = L 9ip / p F and a: 

£ = fL F (* i 1 ] (5) 

Rewriting Eq. (|4]) in the reference frame {e r , eg, e v }, and recalling that the current is the 
density multiplied by the velocity we obtain the well known LL equation: 



-L F |e r X ^H eff + -y^P^j + °&r X (zr X (fieff + ~^ Pj ) } " ( 6 ) 



rfe r , f _ / ^» fcT \ _ / _ ( ^ kT , 

The equivalence between the LL equation and the phenomenological Gilbert equation 
gives the relation between the coefficients a and Lp on the one hand, and the Gilbert 
coefficients r\ and the gyromagnetic ratio 7 on the other hand: 

a = rf{M s ^ 

Lf = M s (l+a 2 ) 

Inserting Eq. ^ and Eq. Eq. Q into the continuity equation Eq. ([I]) leads to the 
Fokker-Planck equation in the configuration space T, x : 

j f 

= V s . jp F e r x (L F # e// + £V s p F ) + ap F e r x (e r x (l f # £// + W s /)) } . (8) 

where D = LpkT/p F is the diffusion coefficient. This result is well known [TTJUH]. The goal 
is to generalize the description to non-uniform ferromagnets (i.e. performing the extension 
of the configuration space to the neighbors T, x ±$ x , in order to describe transfer of magnetic 
moments (beyond Eq. ([!])) and transfer of energy. The objective of the next section is to 
simplify the problem by defining fist the two channel model for magnetic excitations. 



B. The two spin-channel model: reduction of the ferromagnetic degrees of freedom 
to a bivaluted variable. 

The concept of "spin-channels" is inspired from spin-dependent transport studies (or 
"spintronics" ) , for which the electronic spin is a bivaluated degree of freedom with s = ±h/2 
corresponding to the two states "| and 4- for fixed quantization axis. 

We will show in the following the a ferromagnetic bivaluated variable can be defined by 
reducing the continuous degrees of freedom of the magnetization to a bivaluated variable, 
i.e. to an effective spin. This reduction is performed by the integration of the relevant 
quantities over the two hemispheres E| and E^, such that E x = E| U Ej; (as shown in Fig. 
1). We can then define the chemical potentials of the two hemispheres by: 

pL^(x) — / p(#, if; x)dQ, (9) 
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The total chemical potential is given by the sum over the two hemisphere /i F = /j, f + /if , 
and the difference defines a "pumping force" [22] A/x = fi F — /xf. The number of magnetic 
moments n F per unit length Sx for each hemisphere reads: 

n F {x) = [ p F ((p,6;x)dtt, (10) 

so that the density of magnetic moments per unit length is n F = n F + n F and the difference 
between the hemispheres is An F = n F — nf. 
Let us define the scalar J F by the expression: 

JP( X ) = [ divxJ F dQ, (11) 

and the difference 



AJ F (x)= ! divnJ^dn- f div^P dtt. (12) 
From Eq. ([!]), Eq. (10) and Eq. (11) we have: 



^ _ <tF _ 



dt 2 
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where the last equality in the left hand side is obtained by noting that the total numbers 
of magnetic moments on the sphere no = n^ + n^ is constant, so that dn F /dt = J F + J F = 0. 

Note that AJ F is the flux of the magnetic moments flowing from the hemisphere S| to 
the other hemisphere Sj, and the "pumping force" A/x is thermodynamically conjugate to 
the "flux" ^f- in the sense that the product V^^ — AJ F A{i/1 is the ferromagnetic power 
exchanged between the two hemispheres. 

In the same way as above, the ferromagnetic energy u F (8, tp; x) is integrated over the two 
hemispheres. The reduced variable reads: 



and from Eq. §2§ we have: 



U^(x) = I u F ((p,9;x)dQ, (14) 



fr J " +n/2 d5) 
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where 11/2 = f s rdVL = J E ^ rdQ, and the energy currents are defined on each point 
x by the relation: 

j"(x) = [jiv^dn. (16) 



If the system is adiabatic total energy Uq =U^ + Ui on the sphere E x would constant and 
we would have, as for the current of magnetic moments: 



' (17) 



't 2 



U _ AJ U 
I - 2 ' 



where AJ U = J}* — J^ 1 . In conclusion, even for an isolated and uniform ferromagnetic 
particle (i.e. adiabatic ferromagnetic system), a current of energy is flowing from one hemi- 
sphere to the other at any point x. In other terms, an effective "spin-flip relaxation" has 
been defined for ferromagnetic excitations. 



C. Transport equations for heat current along the x direction. 

A ferromagnetic wire can be modeled by a succession of segments of thickness 5x and 
section unity [6] . At each position x, the magnetization of volume ±5x is described in the 
configuration space S^, with the corresponding chemical potential. The configuration space 
of the total ferromagnetic wire of length / is then the continuous limit (5x — > 0, iV — > oo and 
N5x = I) of a chain of N ferromagnetic configuration spaces ... U T, x sx U T, x U T, x+ s x U ... 
(see Fig. 2). 

Let us first consider an open system that is able to exchange magnetic moments with its 
environment [Bj. In this case, a mechanism of transport of magnetic moments would take 
place in the x direction, and the conservation equations for magnetic moments through the 



space x (Eq. (13)) would take the same form as in the case of spin-dependent transport in 
the two channel model for electric conductors: 

drif _ djf • p 



dt dx 



dt dx 

— - — sr + V 



(18) 



where ip F = AJ~ F / (25 x) is the effective spin-flip relaxation, in the sense that the flux ip F is 
the velocity of the transformation of "spin up" \ into a "spin down" \, under the action of 
the chemical affinity A/i F 0HSJ ESI HO 123 I 
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FIG. 2: Illustration of the configuration space £ for a non- uniform ferromagnet of length I = N5x 
and section unity, defined as a chain of N spheres: S = ... U H x sx UEjU T, x+ $ x U ... at the 
continuous limite (N — > oo, Sx — >• 0). 

However, in a ferromagnet, magnetic moments are localized so that a spin- injection mech- 
anism carried by the magnetic excitations is not expected. The magnetic moments cannot be 
transmitted from the configuration space E| to its neighbors E^ ±1 and the local configuration 
space is closed from the point of view of the ferromagnetic degrees of freedom. 

In contrast, since the system under interest is not adiabatic energy is transmitted from 
the configuration space E| to its neighbors E| ±1 . Typically, this situation corresponds to 
ferromagnetic resonance or heating performed at an extremity of the ferromagnetic device 
[SHU]. Space-dependent and spin-dependent heat currents J^(x) are then produced and 



flow throughout the sample and through the interfaces. Furthermor, as shown in Eq. (15), 
the energy flux AJ~ U is produced inside the configuration space E x , flowing from the subspace 



dUf __ 




dt 


dx 


dU l _ 




dt 


dx 



E| to the subspace E^. Generalizing Eq. (15) with the contribution of the neighbors Ti x ±i, 
we have: 

' - v lh - K/2 

(19) 

+ ip th + K/2, 

where ip th = is due to the relaxation between the two channels. The heat current 
generated by the energy carriers is given by the relation: 

3% = 3% ~ n^x (2°) 

We are here interested in the case = 0. The diffusion equations for each channel can 
then be deduced [7j: 
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D. Conclusion 

A non-equilibrium thermodynamic approach has been used in order to establish a two 
channel model applied to macroscopic ferromagnetic excitations. The reduction of the ferro- 
magnetic degrees of freedom from two continuous coordinates {9, <p} to a bivaluated scalar 
variable £ (and a quantification axis) is explicitly performed, and applied to the relevant 
physical observables, i.e. the density of magnetic moments, the energy density, and the 
corresponding currents. The complexity of the transport properties of magnetic moments 
and magnetic energy (including diffusion) can be reduced to the kinetic equations of two 
states f and 1 at each position x in space. This method allows the spintronics concept 
of spin-injection to be generalized to any kind of magnetic excitations. In particular, it is 
shown that thermal spin-injection can be performed without magnetic carriers, but with 
spin-dependent heat currents only. The validity of the usual technics of spintronics has been 
extended to the case of ferromagnetic excitations (quantal or classical) occurring in electric 
conductors or electric insulators. 
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